D and Dg meson spectroscopy 
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Results are presented for the low-lying spectrum of D and Ds mesons calculated in lattice QCD 
using 2-1-1 flavor Clover- Wilson configurations made available by the PACS-CS collaboration. For 
the heavy quark, the Fermilab method is employed. The main focus is S- and P-wave states of 
charmed and charmed-strange mesons, where previous lattice QCD results have been mostly from 
quenched calculations. In addition to the ground states, some excited states are extracted. To check 
the method, calculations of the charmonium spectrum are also carried out. For charmonium, the 
low-lying spectrum agrees favorably with experiment. For heavy-strange and heavy-light systems 
substantial differences in comparison to experiment values remain for some states. 

PACS numbers: ll.15.Ha, 12.38.Gc 



I. INTRODUCTION 

The spectrum of charmed-strange mesons^ contains a 
number of well-established states notably the "S- 
wave" states of quantum numbers 0~ (where J is the 
spin and P is parity) and I~, the Dg and the D* and the 
"P-wavc" states with quantum numbers 0+ (D*o(2317)), 
1+ (i:»,i(2460) and D,i(2536)) and 2+ (£1:2(2573)). In 
the heavy quark limit, these states would form three 
mass-degenerate multiplets and are characterized by the 
total angular momentum j of the light quark. Within 
many quark models (see for example 0), the physical 
states corresponding to the doublet with = i, the 
Z?*o(2317) and Dsi{2460) were expected to have consid- 
erable width and masses above the DK and D*K thresh- 
olds. In experiment, both states are below threshold and 
narrow. 

In addition there are a number of Dg states which have 
been observed more recently. These are the £'*,(2700) 
(observed by both BaBar [3 and Belle jsj), the 
£1:^(2860) (observed by BaBar [1,||), the Dgj{3040)) 
(observed by BaBar Q) and an unconfirmed state pre- 
viously observed by SELEX Q, the Dsj{2632). While 
the D*i{2700) is commonly believed to have quantum 
numbers = 1^, there are several possibilities for the 
other states which are not ruled out by experiment. The 
ZJ^j (2860) has natural parity and is most often identified 
with a 3~ state 041], while some still argue the possibil- 
ity of a 0+ identification The (3040) has unnat- 
ural parity and is commonly interpreted as either a 1+ 
0-[ilni or a 2- state 0. 

Until recently, the number of observed charmed meson 
states was more limited. Again the "S-wave" states of 
quantum numbers Q~ {D) and 1~ (£)*) as well as the 
"P-wave" states with quantum numbers 0+ (150(2400)), 



^ Following Particle Data Group usage [l[ Ds mesons are referred 
to as charmed-strange and only D mesons as charmed. 



1+ (L>i(2420) and £'1(2430)) and 2+ (£)^(2460)) are well- 
established In addition, the Particle Data Group [l| 
lists the £'*(2640) seen in Z-decays [l2| which lacks con- 
firmation. More recently BaBar observed several new 
charmed mesons [l3l | whose quantum numbers are a mat- 
ter of active discussion [1Jt13 • 

Lattice QCD (LQCD) provides the possibility to elu- 
cidate the spectrum without resorting to model assump- 
tions. To reach this goal, several systematic sources of 
uncertainty have to be controlled. Of particular impor- 
tance is that full dynamical simulations are necessary 
with up and down quarks light enough so that an ex- 
trapolation to the physical point is not needed or can 
be done in a meaningful way. Also, there should be an 
extrapolation to the continuum limit. Such calculations 
have been achieved for the light-quark ground state me- 
son and baryon spectrum (see, for example, fioj). 

Dealing with heavy charm or bottom quarks on the lat- 
tice introduces complications of its own, and so far a sim- 
ilar precision for the spectrum of charmed and charmed- 
strange states has not been attained. Most previous lat- 
tice calculations have been within the quenched approxi- 
mation or at unphysically large pion masses. Advances in 
lattice simulations enable us to revisit this problem, re- 
moving one deficiency of previous calculations by using 
dynamical gauge ensembles where pion masses are close 
to the physical pion mass. 

Among the low-lying D and Dg mesons accessi- 
ble to current lattice calculations, the £':o(2317) and 
-Dsi(2460), whose properties can not be explained within 
most quark models, are of particular interest. Possible 
explanations for their unanticipated mass and width have 
been reviewed in [13, [Hj. In early LQCD calculations, 
the ground states in the 0+ and 1+ channels were often 
found to be quite a bit heavier than the experimental res- 
onances (for a short summary of previous results please 
refer to Section ITV]) . This has further sparked specula- 
tions about the nature of these states and provides mo- 
tivation for this study. 

The next section provides some details of our calcula- 
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tional setup. There is also a short description of the pro- 
cedure used to tune the charm quark mass. To test our 
setup, the spectrum of charmonium states below multi- 
particle threshold was calculated and is discussed in Sec- 
tion mil Section IIVI presents our results and discussion 
of the charmed and charmed-strange mesons. Finally, in 
the last section, a summary of our findings is presented. 
Also included is a short description of the role of scat- 
tering thresholds for the calculation. Some preliminary 
results have already been presented in [22[ . 



II. CALCULATIONAL SETUP 

A subset of the dynamical 24-1 fiavor configurations 
generated with the Clover- Wilson fermion action by the 
PACS-CS collaboration [2^, spanning sea-quark pion 
masses from 702 MeV down to 156 MeV, are used in 
this work^. The number of lattice points is 32^ x 64 
for all ensembles, and using the sing le lattice spacing of 
0.0907(13) fm, as determined in [23 this corresponds to 
a box size of roughly 2.9fm in spatial direction and 5.8fm 
in time direction. Table |I] shows the parameters for the 
runs. With exception of the ensemble #5 we use ~ 200 
configurations on each ensemble. 



Ensemble 




l^u 1 d 




#configs DjDs 


1 


1.52617 


0.13700 


0.13640 


200/200 


2 


1.52493 


0.13727 


0.13640 


-/200 


3 


1.52381 


0.13754 


0.13640 


200/200 


4 


1.52327 


0.13754 


0.13660 


-/200 


5 


1.52326 


0.13770 


0.13640 


200/348 


6 


1.52264 


0.13781 


0.13640 


198/198 



TABLE I: Run parameters for the PACS-CS lattices [13]. The 
number of lattice points is 32"^ x 64 for all ensembles. All gauge 
configurations have been generated with the inverse gauge 
coupling /3 — 1.90 and the light quark clover coefficient Csm — 
1.715. The number of configurations indicated in the table 
corresponds to the number used rather than the total number 
available. The quantity cim denotes the heavy quark clover 
coefficient. Please refer to Section Til Dl for the determination 
of cit' . 



A. Source construction 

The low-lyin g sp ectrum is extracted using the varia- 
tional method j27l [28j . For each combination of quan- 



^ After work on this project started dynamical simulations with 
even smaller pion masses, for the first time at [24j and around 
l25l |26(I the physical point, were reported. 



tum numbers (where J is the Spin and P is parity) a 
matrix C{t)ij of interpolators is constructed 

C(%=5]e-*^"(0|O,|n)(n|O]|0). (1) 

n 

Subsequently, the the generalized eigenvalue problem is 
solved for each time slice 

Cit)^'^''^ = X^^\t)C{t„)ib'^''\ (2) 

AW(i)occ-*^^'(l -1-0(0-*^^'=)). 

At large time separation only a single mass contributes 
to each eigenvalue. AEk is in general given by the en- 
ergy difference between the energy level in consideration 
and the neighboring level. For more details the reader is 
referred to the discussion in (29j . For the determination 
of the charm quark hopping parameter we will use pro- 
jections to the lowest few lattice momenta, while masses 
are obtained from a projection to vanishing spatial mo- 
mentum. In addition to the eigenvalues, the eigenvectors 
provide useful information and can serve as a fingerprint 
for a given state. 

The basis is constructed from two types of sources. The 
first type arc simple Jacobi-smeared [s^, HH Gaussian- 
shape sources = (S q)x 

N 

S = M So with M = ^ 

3 

H{ri, m)^Y^ (Uj {n, 5 [n + j, rri^ (3) 

Here, H is the spatial hopping term containing the link 
variables Uj (n, t) and the smearing has two parameters 
K and N. In addition, we also use derivative sources Wd^ 

Wd, ^D,S , (4) 

D.i{x,y) = Ui{x,0)S{x + i,y) - Ui{x - 5[x - i,y) . 

The parameters for the Gaussian-shape sources arc k = 
0.22 and N = 25. For the derivative sources the num- 
ber of iterations is slightly larger A^ = 30 while using 
the same k. The full list of interpolating fields used in 
the simulation can be found in Appendix lA II Similar 
source constructions have been used for both quarko- 
nium [H, nil and light meson [33 - l36j spectroscopy. To 
minimize the calculation of expensive propagators, we 
use both Gaussian and derivative sources for the charm 
quark but only Gaussian sources for the strange and light 
quarks. This has an additional implication: The inter- 
polators with derivative sources are constructed from a 
derivative which is not symmetrically applied to both 
quarks and these interpolators do not have good charge 
conjugation quantum number. This is relevant in the 
case of charmonium, where the same basis is used. For 
further comments regarding charge conjugation, please 
refer to Section Hill 
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B. Quark propagators 

For the construction described in the previous section 
we need a single full propagator for the strange (light) 
quark and 4 full propagators (one Gaussian and three 
derivatives) for the charm quark. These propagators are 
calculated for 8 source time slices on each gauge con- 
figuration. The source locations are chosen randomly 
within the time slice. For the calculation of the light and 
strange quark propagators the dfl_sap_gcr inverter from 
Liischer's DD-HMC package [13, HI] is employed. For 
the charm quark propagators the corresponding inverter 
without deflation is used. The large number of sources 
needed for this work makes the use of a deflation inverter 
extremely efficient, particularly for light quarks on En- 
semble 6, with a pion mass of 156MeV. 



C. Fitting methodology 

This section presents the methodology for fits to both 
single correlators as well as eigenvalues from the varia- 
tional method. All of these fits take into account au- 
tocorrelation in Euclidean time t. The estimate of the 
covariance matrix Cov(i, t') employs a single-elimination 
jackknife 



Cov{t,t') 



N ■ 



N 



1 ^ 

i^(A«(t)-A(t)) (a«(0 
1=1 



lit'] 



(5) 



where N denotes the number of configurations and the 
bars denote averages. To obtain a stable estimate, the 
window for the inversion of the covariance matrix is re- 
stricted to the chosen fit range. The inverse covariance 
matrix is estimated once on the ensemble average and 
used for each jackknife block. This method has been 
dubbed "jackknife reuse" in [s^. To determine the win- 
dow in Euclidean time for performing fits we calculate 
effective masses 



aA'Icff = In 



A(')(t) 
A(^)(t+1) 



(6) 



and eigenvector components of the regular eigenvalue 
problem 

C(to)-^C(t)C(to)-^>^" = AW(t)V^W', (7) 

and choose as our fit range the interval where both the 
effective masses and the eigenvector components are con- 
stant within their statistical uncertainties. Then two pa- 
rameter fits are carried out within this window. As a 
cross check, 4 parameter fits with two exponentials and 
a larger fit range are performed. In some cases, the 4 pa- 
rameter fits are unstable. Where stable, the 4 parameter 
fits lead to results that are consistent within errors with 
the strategy described above. The dependence of the 



results on the time slice was investigated and, where 
necessary, the reference time to was increased until the re- 
sults were qualitatively unaffected by a further increase. 
Appendix I A 21 collects various fit results, fit intervals and 
the associated x^/d.o.f.. 

Due to the method employed for the heavy quark, the 
energy levels extracted should not be identified with the 
physical energy levels directly. Instead, as final results 
we always quote energy differences with respect to the 
spin- averaged IS states, given by 



M- 



{Md^+3Md*J/4, 

MTy= {Md+3Md*)/4, 
A%= (M^^ +37\/j/*)/4. 



(8) 



for the case of heavy-strange mesons, heavy-light mesons 
and charmonium respectively. As the measurements are 
on the same gauge ensemble and therefore correlated, 
taking into account this correlation leads to reduced sta- 
tistical errors for the mass differences. 



D. Heavy quark action and quark mass tuning 

To determine the mass parameter for the heavy charm 
quark we use the Fermilab method j40| in the form em- 

ted by the Fermilab Lattice and MILC collaborations 
for their efforts involving quarkonium 42 1. This sec- 
tion describes this method in brief and points out where 
our method differs from [4l|. In addition to this sim- 
plest prescri ptio n, an improved Fermilab action has been 
constructed [43| and preliminary results are encouraging 
(4^ . Within this approach, the charm quark hopping pa- 
rameter Kc is tuned to the value where the spin-averaged 
kinetic mass [Mo^ -f SM^i* )/4 assumes its physical value. 
In this simplest formulation the heavy quark hopping pa- 
rameter ce = cb = citi* is set to its tadpole improved [45| 
value -ij , where uq is the average link. For simplicity, uq 

is estimated to be the fourth root of the average plaque- 
tte, unlike in [4lj, where it is given by the Landau link. 
The lattice dispersion relation takes the general form [4l| 



E{p) = Ml 



2A/2 



(9) 



where the momentum p = and L is the spatial extent 
of the lattice. To determine the kinetic mass M2 the six 
lowest lattice momenta are used and averaged over all 
possible vectors n. These are denoted symbolically as 
000, 100, 110, 111, 200 and 210. Larger momenta are very 
noisy and do not help to constrain the fits further. As 
there are not enough points to constrain a fit to the four- 
parameter form of Equation [SI two simplified methods 
are considered: 

1. neglect the term with coefficient W4 and fit Mi, 
M2 and M4. 



2. fit E'^iji) and neglect the (p^)^ term arising from 
the mismatch of Mi, AI2 and M4 



Ml 
M2 



Mia^W4 



(10) 



in physical units has two contributions and the last row 
of Table HI] lists both the statistical uncertainty and the 
uncertainty from setting the lattice scale. This second 
error from the conversion to physical units is dominant 
for method 1. 



For both methods, two values of the charm quark hop- 
ping parameter close to the physical value are used and 
the results are interpolated linearly. To determine the 
energy levels, a 2 x 2 matrix of interpolating fields cor- 
responding to interpolators A and B in Appendix lA II is 
used in both the = 0~ and 1~ sectors. We choose 
to = 3 for these fits. Table HI] shows the results for 
the spin-averaged ground state from the first method. 
For illustration an example fit from this method with 
Kc = 0.128 is shown in Figure [H 





Kc = 0.128 


Kc = 0.127 


Ml 


0.8633(5) 


0.8931(5) 


Ma 


0.9337(73) 


0.9716(76) 


M4 


0.8638(274) 


0.8855(284) 


M2 

Ml 


1.0815(86) 


1.0878(88) 


M2[GeV] 


2.0315(158)(291) 


2.1137(166)(303) 



TABLE II: Fit parameters obtained with our first method for 
two different values of Kc- The values in the last row are in 
GeV, while all other values are in lattice units. The first error 
on the kinetic mass M2 is statistical while the second error is 
from the scale setting. The results for M4 are not used in our 
setup. 



0.98 p 
0.96- 
0.94- 
0.92- 
0.9- 
0.88- 
0.86- 



o Lattice data 

— fit result - method 1 



FIG. 1: Example fit using method 1 with Kc — 0.128 

The three parameter fits are stable and both the rest 
mass Ml and the kinetic mass M2 are well-determined. 
The mass M4 is not used in our analysis but it is en- 
couraging that the result is consistent with the expec- 
tation Ml ^ M2 ^ M4. The error on the kinetic mass 





Kc = 0.128 


Kc = 0.127 


Ml 


0.8634(5) 


0.8932(5) 


Ml 


1.0889(116) 


1.0955(118) 


M2[GeV] 


2.0454(215) (293) 


2.1293(227)(305) 



TABLE III: Fit parameters obtained with our second method 
for two different values of Kc. The values in the last row are 
in GeV, while all other values are in lattice units. The first 
error on the kinetic mass M2 is statistical while the second 
error is from the scale setting. 




0.127 



FIG. 2: Linear interpolation to determine the charm quark 
hopping parameter Kc- The horizontal line corresponds to the 
physical value for {Mds +3Afoj)/4. The red curve (squares) 
is a shift taking into account the slightly unphysical strange 
quark mass. 

The results from method two can be found in Table 
mil This second approach leads to slightly larger values 
for the kinetic mass M2. The results of the two meth- 
ods are however consistent within the statistical uncer- 
tainties, which are slightly larger for the second method. 
With both methods M2 is determined sufficiently well, in 
the sense that the error due to the lattice scale setting is 
dominant. The calculation of the heavy-light and char- 
monium spectra in Sections IIIII and IIVI uses the estimate 
of Kc from fits with the second method. For this, a plot of 
the results for the spin-averaged ground state is shown in 
Figure [2] To account for the slightly unphysical strange 
quark mass used in the simulation we also calculated the 
spin-averaged rest mass with a partially quenched strange 
quark mass corresponding to the estimate of the physical 
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strange quark mass from [23| • The red curve (squares) in 
Figure [5] results from shifting our original results to ac- 
count for this difference. One can then proceed to read 
off the value Kc = 0.12752 with which the simulations to 
determine the low-lying spectrum are performed. 

Having tuned the charm quark hopping parameter us- 
ing Ds mesons, a stringent test of our setup can be per- 
formed by calculating the low-lying charmonium spec- 
trum, where there are several well-established states be- 
low multiparticle-threshold. 



III. CHARMONIUM RESULTS 

In recent years many new charmonium resonances have 
been observed, primarily at the B-factories. While the 
agreement between the well-established low-lying states 
and quark model calculations (see for example 0]) is 
quite good, many of the newly observed X, Y and Z 
states do not fit within a simple quark model picture. For 
a recent review regarding these charmonium-like states 
please refer to [461 . For an emphasis on recent results 
see j47l - l49| . The charmonium states considered in this 
section are all well established and considered to be pre- 
dominantly regular quark-antiquark states with masses 
well below the DD and DD* thresholds. This makes the 
calculation of the low-lying charmonium spectrum a good 
test of the heavy-quark methodology. Where appropri- 
ate, we somewhat sloppily refer to those states with their 
spectroscopic nomenclature and call them S- and P-wave 
states respectively. 

As mentioned in Section III Al some of our interpola- 
tors, which were chosen with heavy-light mesons in mind, 
mix different charge conjugations. In the case of the rjc 
(J^c = 0-+), the J/^ (J^^ = 1 — ) and the Xco chan- 
nels {J^'~^ — 0'^'^) the corresponding admixtures would 
come from exotic states, which are expected to play a role 
higher in the spectrum but not in the vicinity of the low- 
lying states extracted in this section. In the case of the 
1+ channel(s) the basis is restricted to those interpola- 
tors with good charge conjugation and for the interpola- 
tors containing the Xc2 we just extract the ground state, 
which again should be safe from possible contamination. 

Figure [3] shows the chiral behavior of the charmonium 
hyperfine-splitting Mj/q, — M^i^ compared to the experi- 
mental value. The errors in the plot are purely statisti- 
cal and errors from the tuning of the charm quark mass 
are non-negligible'^. Within our approach, the hyperfine- 
splitting is expected to show significant discretization er- 
rors [44I which we can not further investigate, as the 
library of PACS-CS lattices used contains only one lat- 
tice spacing. This result indicates that discretization ef- 
fects of a similar magnitude may be expected in other 
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^ Please refer to Section IIVI for an estimate of this uncertainty in 
the case of the Ds hyperfine-splitting. 



FIG. 3: IS hyperfine-splitting for charmonium compared to 
the physical splitting. The error is statistical only and the 
leading discrepancy is caused by the lack of a continuum ex- 
trapolation. For further discussion, please refer to the text. 



observables. In addition, disconnected diagrams which 
are neglected in our simulation will have a small effect 
which may be non-negligible for the hyperfine-splittings. 
In [soj these contributions to the charmonium hyperfine- 
splitting have been calculated for a reasonable fitting 
Ansatz. They are found to be negative and therefore 
can not explain the difference of our results compared 
to experimental values. Our results also show a small 
but statistically significant dependence on the strange 
sea quark mass. Notice however that our charm quark 
tuning was performed for the ensemble with lightest sea 
quark mass and = 0.1364. For all other charmonium 
mass differences the dependence on the strange sea quark 
mass is very mild and there is no statistically significant 
difference between the two ensembles at Ku/d = 0.13754. 

Next consider the spin-triplet P-wave states, which are 
the ground states with quantum numbers 0^^, 1+^ and 
2"^+ corresponding to the XcOj Xci and Xc2- Figure 2] 
shows the difference between the masses of the individ- 
ual states and the spin-averaged ground state mass for all 
ensembles. At the lowest quark mass the ground state P- 
waves calculated at finite lattice spacing show good qual- 
itative agreement with experiment, with the XcO ground 
state showing the largest discrepancy. While our results 
at the three heaviest and at the lowest quark masses sug- 
gest that the chiral extrapolation is within errors well de- 
scribed by the leading order term linear in the pion mass 
squared, the data from the ensemble at Ku/d = 0.13770 
deviates from this behavior. This peculiarity is found in 
most of the observables we investigate. Note that a sim- 
ilar behavior can be seen in [5l| where the first positive 
parity excitation of the Nucleon, the Roper resonance is 
calculated on the same gauge configurations. In [5l| this 
is interpreted as "significant chiral curvature" . Our data 
arc not suggestive of this explanation. Without addi- 
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FIG. 4: Chiral behavior for the single charmonium P-wave 
states compared to experiment. As for all plots, the numbers 
are mass differences with respect to the spin-averaged ground 
state. The error displayed is statistical only. The uncertainty 
in the lattice scale implies that all points could be moved 
uniformly by ~1.4%. 
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FIG. 5: Splitting Mj-p — Mjg between the spin-averaged IP 
triplet state and the spin-averaged ground state. The error 
displayed is statistical only. The uncertainty in the lattice 
scale implies that all points could be moved uniformly by 
R^l.4%. 



tional simulation data we refrain from interpreting this 
behavior and will instead quote as our final results those 
at the lightest pion mass m.^- w 156McV. The effect of the 
remaining small chiral extrapolation would be negligible 
compared to current statistical and other systematical 
uncertainties. 

In addition, we also form the combinations (see [4l| 
and references therein) 



IP 



M, 



Spin— Orbit 



(5M^ 



3Mv 



2Mv 



2NK 



(11) 
(12) 

(13) 



il/j-p is the Spin- averaged mass of the P wave triplet. 
Mspin-Orbit IS Sensitive to the spin-orbit interaction 
while Mxensor IS Sensitive to the spin-spin interaction 
j4l| . Figure [5] shows the splitting between the spin- 
averaged ground state mass Mj^ and Mjp. While we 
are unable to perform a continuum extrapolation and 
small discrepancies to experiment are expected, the fig- 
ure suggests that the spin-averaged quantities are well 
reproduced within our approach. 

The results for Mspin-OrUt and Mxensor are plotted 
in Figure ini The upper panel shows the Spin-Orbit split- 
ting of the IP states, which is slightly underestimated 
in our simulation, similar to the hyperfine-splitting. The 
lower panel shows the tensor splitting of the IP states. 
Compared to experiment we obtain the correct sign and 
a similar magnitude with the value off about 30%. 

So far the discussion was restricted to properties of 
ground states. In the 1 and channels our basis 
allows for a determination of excited energy levels. While 
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FIG. 6: Spin-orbit splitting (upper panel) and tensor splitting 
(lower panel) for the triplet of IP states. The error does 
not inchide the error due to the tuning of the charm-quark 
hopping parameter Kc and is statistical only. The uncertainty 
in the lattice scale implies that all points could be moved 
uniformly by ~1.4%. 



0.8 
0.7 
> 0.6 

o 

^0.5 
o 

I 0.4 
^ 0.3 
I 0.2 
0.1 



O lattice K_,= 0.1364 
A lattice K =0.1366 

s 

O IS - 2S splitting experiment 



?0^ 



0.1 



0.2 



0.3 



0.4 



0.5 



[GeV'] 



FIG. 7: Splitting - Mjg between the 2S and IS spin- 
averaged states. The error displayed is statistical only. The 
uncertainty in the lattice scale implies that all points could 
be moved uniformly by ~1.4%. 



0.16r 



0.14- 



o lattice K^= 0.1364 

0.12- A lattice K^= 0.1366 

o 2S hyperfine experiment 

0.1- 

0.08- 
0.06- 
0.04- 



0.02- 



-0.1 



0.1 



0.2 0.3 
M ^ [GeV^] 



0.4 



0.5 



these states are more noisy than the ground states and 
the fit ranges (see tables in the Appendix lA 2[) are more 
limited, at least the 2S states are well determined from 
our charmonium data. Figure [7] shows the splitting be- 
tween the 2S and IS spin-averaged states. Again the re- 
sult agrees qualitatively with experiment suggesting that 
our approach is well suited to reproduce spin-averaged 
quantities. 

Figure [5] shows the 2S hypcrfine-splitting. Unlike the 
IS hypcrfine-splitting this quantity is not determined 
very accurately from our data. Just like in experi- 
ment the 2S hyperfine-splitting is smaller than the IS 
hyperfine-splitting. This again gives us confidence that 
we capture the important physics well. This concludes 
the discussion of charmonium results. Final numbers 
from the ensemble with lightest sea quarks can be found 
in Section |V] along with a table that summarizes the re- 



FIG. 8: Hyperfine-splitting for the 2S charmonium states 
compared to the physical splitting. The error is statistical 
only. 



suits and their respective uncertainties. A table with 
charmonium energy levels can be found in Appendix lA 21 



IV. HEAVY-LIGHT MESONS 

In the heavy quark limit, the spectrum of S- and P- 
wave heavy-light mesons can be classified by the to- 
tal angular momentum j of the light degrees of free- 
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dom . This classification contains one S-wave doublet 
with = i containing the D and D* mesons in case 
of the heavy-light states and the Dg and D* mesons 
in case of the heavy-strange states. The P-wave states 
fall into two multiplets, one with = i"*" and one 

with ~ 1^. Taking a look at the experimental sit- 
uation one would identify the charmed-strange mesons 
-D*g(2317) and Dsi{24:60) as the charmed-strange ver- 
sion of the P-wave heavy quark doublet with j ~ and 
the Z?si(2536) and the Z)*2(2573) as the charmed-strange 
version of the multiplet with J = §• While the mass of 

states 



states corresponding to j' 



and J — 



in the heavy mass limit is quite close to the values ex- 
pected from the quark model (H , the states corresponding 
to the heavy-strange = doublet are substantially 
lighter [J , putting them below the D + K and D* + K 
thresholds. Even more puzzling however is a comparison 
to the spectrum of heavy- light D mesons, where similar 
multiplets are expected in the heavy quark limit. In par- 
ticular, the charmed-light version of the = i"*" doublet 
is given by the D5(2400) and the £'i(2420). This corre- 
sponds to a mass-splitting between these states and the 
spin-averaged S-wave states which is much larger for D 
mesons than for Dg mesons. Such a behavior is not ex- 
pected in simple potential models. Furthermore, chiral 
perturbation theory for static-light mesons suggests the 
opposite (5^ . 

Before proceeding to report results, let us discuss 
previous lattice calculations of charmed and charmed- 
strange mesons focusing in particular on results for the 
~ doublet. Lewis and Woloshyn [s^ calculated 
the spectrum of S- and P-wave heavy-light mesons in lat- 
tice NRQCD using quenched gauge configurations. At 
the time of this study, the = 0+ states were not 
yet measured. While their results for the heavier = | 
doublet agree within errors with the experimental results, 
the Dq, and the = \ Dsi were predicted at masses 
much larger than the experimental masses. Another early 
calculation using lattice NRQCD was performed by Hein 
et al. [H^l- In [H^, [HB] the UKQCD collaboration reports 
on results using a relativistic clover-type action based on 
the Fermilab method |4^. Results from the studies (53 - 
[sgI I along with a calculation in the static limit on = 2 
dynamic gauge configurations have been presented by 
Bali (s^]. All these results share the common feature 
that the ground state in the = 0"*" channel is esti- 
mated at a mass exceeding the mass of the experimental 
D*o(2317) by roughly 130-200 MeV. Furthermore, where 
calculated, the splitting between the S-wave states and 
the 0"'' P-wave states for heavy-strange mesons was deter- 
mined to be larger or of the same size than for heavy-light 



* Note that we use small j when referring to the heavy quark limit 
and J when referring to the Spin away from the heavy quark 
limit. 



0.16 



0.15 

> 

^0.14 
a 

9 0.13 



0.12 



0.11- 



o D hyperfine splitting K^= 0.1364 
O experiment 



0.16 



0.1 



0.2 0.3 
M ^ [GeV^] 



0.5 



0.15- 

> 

^0.14 
a 

a 0.13 



0.12 



0.11 



O hyperfine splitting K^= 0.1364 
A hyperfine splitting K^= 0.1366 
O experiment 



0.1 



0.2 
M 



0.3 



0.4 



0.5 



2 7 

. [GeV] 



FIG. 9: Hyperfine-splitting for heavy-light (upper panel) and 
heavy-strange (lower panel) mesons compared to the physical 
splitting. The errors are statistical only and the leading dis- 
crepancy is caused by the lack of a continuum extrapolation. 
For further discussion, please refer to the text. 



states. 

In addition to the above simulations, the UKQCD col- 
laboration also presented results from both quenched and 
one dynamical ensemble in (ssj . While they report a 
somewhat smaller discrepancy with respect to experi- 
ment, their scale setting results in a quite large value 
of 0.55fm for the Sommer scale rg. Typical values for the 
Sommer scale from recent dynamical simulations however 
suggest a significantly smaller value tq « 0.46 . . . 0.49fm. 
In particular the scale setting for our simulations, as de- 
termined by the PACS-CS collaboration, corresponds to 
ro = 0.4921(64)(-h74)(-2) [H. A smaller value of tq 
results in larger mass splittings. The scale setting itself 
is intimately related to the determination of the physical 
point, which again emphasizes the need for light dynam- 
ical quarks. 

More recently, the xQCD collaboration presented pre- 
liminary results from both quenched [s^ and dynamical 
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[6 Of lattices using relativistic quarks and overlap valence 
fermions. Their initial results for both quenched and 
dynamical simulations suggest good agreement between 
lattice results and the low-lying charmed-strange mesons. 

The rest of this section presents results for the spec- 
trum of heavy-strange and heavy-light mesons. The final 
numbers and an overview of the results can be found in 
Section El 

Figure [9] shows the hypcrfine-splittings for the D (up- 
per panel) and Ds (lower panel) mesons. Similar to Char- 
monium, the hyperfine-splitting turns out to be some- 
what to small, which we attribute partly to non-negligible 
discretization effects. This discrepancy is somewhat 
smaller for heavy-light states than for charmonium. Like 
for charmonium the hyperfine-splitting on the ensemble 
with Kg = 0.1366 is slightly smaller than for the other 
ensembles. Another non-negligible source of error is the 
uncertainty in the determination of Kc- From the tuning 
runs, this uncertainty is estimated to be roughly 3.2 MeV 
for Dg mesons and we assume a similar dependence for 
D mesons. In the case of Dg mesons where statistical 
uncertainties are small, neglecting the data point at our 
second lowest quark mass, the remaining data could be 
fit with a good x^/d.o.f. by a fit linear in the pion mass 
squared. Although the results have much larger statisti- 
cal errors, chiral effects seem to be more significant for 
D mesons. 

Figure shows results for the D (left panels) and 
Dg (right panels) P-wave states. As for charmonium, 
mass differences relative to the respective spin-averaged 
S-wave state are plotted. The top row shows the ground 
state in the Dq and D*q [J^ = 0+) channels. In both 
cases, the mass splitting gets smaller as the light sea 
quark masses approach the chiral limit. Unlike suggested 
by experiment, the splitting for the Dq is smaller than 
for the D*gQ. While the discrepancy between our data and 
the 1)^0(2317) is smaller for almost physical light quarks, 
calculation at light sea quarks alone can not resolve the 
apparent puzzle between the experimental results and 
the lattice data. For the D mesons, statistical errors 
are larger. As previously remarked for the charmonium 
hyperfine-splitting in Section IlIIl the data calculated on 
Ensemble 5 (see Table |l| seem to lead to larger mass 
splittings for all channels, while the rest of the data sug- 
gest that the leading order term proportional to the pion 
mass squared dominates the chiral extrapolation. This is 
especially evident for the ZJJ meson. For the final num- 
bers and overview plots in Section [Vj we therefore quote 
the numbers obtained on the ensemble with the lightest 
pion mass. Any sensible functional form for a chiral ex- 
trapolation would alter these results by an amount much 
smaller than the dominant statistical and systematic er- 
rors. 

The middle panels of Figure [TU] show the result for the 
two lowest states in the Di and Dgi {J^ = 1+) channels. 
Both plots are from a matrix of correlators which in- 
cludes interpolating fields corresponding to both positive 
and negative charge conjugation in the mass-degenerate 



case. In [57| the importance of this mixing has been em- 
phasized, but the actual calculation neglected the mix- 
ing. In [53 the mixing has been considered, but only the 
ground state could be resolved. Figure [TT] illustrates the 
effect of this mixing. In particular the standard interpo- 
lators A and D (see Appendix lA 2[) both couple strongly 
to the ground state leading to a much smaller splitting 
than the mixed basis. The results from the 4x4 ba- 
sis in Figure [10] exhibit a mass splitting which is still 
somewhat small compared to the experimental splitting 
but enhanced compared to the analysis without mixing. 
For the Dgi the results for the ground state lead to a 
splitting substantially larger than the £)si(2460), while 
the difference between the first excited state obtained in 
the calculation and the Dsi(2536) is compatible with the 
non-negligible discretization errors one expects. For the 
Di mesons, the ground state displays a somewhat larger 
slope in the chiral extrapolation, while the splitting for 
the first excited state is comparable to the Dgi case. 

Finally the bottom two panels of Figure [10] show the 
ground states of the D2 (left-hand side) and D*2 (right- 
hand side) with = 2^ . For some ensembles we en- 
counter a systematic uncertainty (of a size similar to the 
purely statistical error) with regard to the choice of fit 
range in this channel. On the ensemble with lowest quark 
mass, the DJ2 state is well determined though, and the 
difference from the experiment value is small. The D\ 
ground state is somewhat less well determined and the 
difference from experiment is of a similar magnitude as 
for its D\ partner state in the = I"*" heavy quark 
multiplet. 

For the Dg mesons, it was possible to extract some in- 
formation about excited states in the 0~ and 1~ channels. 
The top panel of Figure [T^] shows the resulting splitting 
between the spin-averaged IS and 2S states. To compare 
with experiment we also plot the single excited Dg and 
D\ states in the middle and bottom panels of Figure [121 
The displayed errors are statistical only. For both 2S 
states, the result for the splitting at the lightest quark 
mass is substantially smaller than on all other ensem- 
bles. On ensemble 4 with Kg ~ 0.1366 the results turn 
out to be somewhat larger than for all other ensembles. 
From experiment, a 1~ excited state, the 13*]^ (2700) is 
known. The result at the lowest quark mass is somewhat 
larger than experiment, but seems reasonable given the 
statistical and systematic uncertainties, which are larger 
for excited states than for the ground states. (Regard- 
ing the importance of possible multihadron thresholds, 
remarks similar to those about the ZJ^q and Dg\ ground 
states apply.) To complete the analysis of the 2S states, 
we also plot the 2S hyperfine-splitting in Figure [13] In 
general one expects that the hyperfine-splitting decreases 
for states higher up in the spectrum. In quark models, 
(see for example [2|) a typical splitting for the 2S states 
would be « 60MeV. For D mesons the BaBar collab- 
oration recently identified candidates [l^ for the corre- 
sponding states and the splitting between those D meson 
states is of comparable size. While the errors are large 
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FIG. 10: Chiral behavior for the D (left panels) and Ds (right panels) P-wave states compared to experiment. The _D(2430) 
(left middle panel) has been slightly offset from the physical pion mass to improve visibility. As for all plots, the numbers 
are mass differences with respect to the respective spin-averaged ground state. The error displayed is statistical only. The 
uncertainty in the lattice scale implies that all points could be moved uniformly by ~1.4%. 



11 



Effective masses diagonal correlators 



1.3r 



1.25- 



f 1.2- 



o interpolator A 
□ interpolator C 
interpolator D 
A interpolator F 



1.15- 



u 1.1- 

*-4— 1 

1 1.05 - 
W 



A ge l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3 6 9 12 15 18 21 24 27 30 



Effective masses from Eigenvalues 



1.3r 

1.25- 



f 1.2- 



1.15- 



^ 1.1- 
^ 1.05- 



1 

0.95 








o ground state 
□ 1st excited state 



X2/d.o.f.= 6.03/6 
X2/d.o.f.= 4.72/6 



12 15 18 21 24 27 30 
t 



FIG. 11: Selected single diagonal correlators (upper panel) 
and lowest two states from a 4 x 4 matrix of correlators (lower 
panel) in the Dsi (J^ = 1^) channel. The data is from En- 
semble 6. The interpolating fields are tabulated in Appendix 



compared to the small splitting, we obtain values consis- 
tent with these expectations. 



V. SUMMARY & DISCUSSION 

We presented results from a dynamical lattice QCD 
calculation of heavy-light mesons. To investigate the ef- 
fects of light dynamical quarks, configurations generated 
by the PACS-CS collaboration [23|, with light quarks cor- 
responding to a pion mass as light as 156 MeV, have been 
used. To test the computational setup, we performed a 
calculation of the low-lying charmonium spectrum, where 
scattering thresholds are less important and previous cal- 
culations within the same framework exist [4l| . The re- 
sults from the ensemble with a pion mass close to the 
physical pion are summarized in Figure 1141 The qualita- 
tive features of the low-lying charmonium spectrum are 
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FIG. 12: Splitting between the IS and 2S Ds mesons. The top 
panel shows the splitting between the spin-averaged IS and 
2S states. The remaining panels show the individual split- 
tings between the excited Ds (middle panel) and D* (bottom 
panel) and the IS spin-average. Only one of these states has 
been clearly identified by experiment. The error displayed is 
statistical only. The uncertainty in the lattice scale implies 
that all points could be moved uniformly by f»1.4%. 
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FIG. 13: Hyperfine-splitting between the Ds 2S states. Errors 
are statistical only. 



well reproduced which gives us confidence in the proce- 
dure used to put the heavy charm quark on the lattice. 
The numerical values for charmonium mass differences 
and estimates of some of the uncertainties are provided 
in Table IIVI With a calculation at only one lattice spac- 
ing it is not possible to make a quantitative estimate of 
discretization effects. However, the expected truncation 
errors for the Fermilab action as presented in Figs. 3 and 
4 of [i^ can provide a rough guide. With the action used 
here, discretization effects in spin-averaged mass differ- 
ences of a few percent could be expected. For more sen- 
sitive quantities, e.g., the hyperfine-splitting, only simu- 
lation can provide a reliable estimate. 
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FIG. 14: Mass differences in MeV for single charmonium 
states compared to experimental values. All splittings are 
with respect to the spin-averaged ground state mass (M^^ + 
3Mj4<)/4. The errors in the plot have been calculated by 
combining the errors listed in Table HVl in quadrature. 

Figure [15] shows the results for heavy-light D and 
heavy-strange Ds mesons as calculated on the ensem- 



TABLE IV: Mass differences in the charmonium spectrum 
in MeV compared to experimental values (calculated from 
Pi). Bars denote spin-averaged values as discussed in Section 
mil For the results of this paper, the first error denotes the 
statistical error and the second error denotes the error from 
setting the lattice scale. In addition there is a non-negligible 
error from the uncertainty in the determination of Kc for all 
spin-dependent quantities. We estimated this error for the 
hyperfine-splitting in the heavy-strange system, but do not 
determine this uncertainty for charmonium. We stress again 
that the gauge ensembles at our disposal do not allow for a 
continuum and infinite volume extrapolation. Consequently 
we expect qualitative but not quantitative agreement. 



ble with the lightest quark masses. In both cases differ- 
ences between the state of interest and the spin-averaged 
ground state are shown. Table IVl provides numerical val- 
ues for these splittings and for the IS and 2S (for the 
Ds) hypcrfinc-splittings. Even at pion masses close to 
the physical pion mass, significant differences in compar- 
ison to experiment remain for the case of Ds mesons. 
This is especially evident for the states corresponding to 
a = doublet in the heavy quark limit. For the 



D mesons, errors are larger and the difference from ex- 
perimental values is more pronounced for the = |^ 

states. The results for the = doublet on the other 
hand are qualitatively different from the heavy-strange 
case. 

While light dynamical quarks improve the agree- 
ment with experiment for heavy-strange and heavy-light 
mesons, considerable discrepancies remain in the case of 
the 0^ and 1+ channels. There are a number of possi- 
ble reasons for this. For one, future investigations will 
have to include a full continuum extrapolation. At our 
lattice spacing, discretization errors are non-negligible. 
However, in the light of previous results for the char- 
monium spectrum, it seems unlikely that discretization 
errors alone are to blame for the observed behavior. An- 
other possible reason for discrepancies is the lack of an 
extrapolation to the infinite volume limit. While the 
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FIG. 15: Mass differences in MeV for D and Ds meson states 
compared to experimental values. The difference is always be- 
tween the state and the spin-averaged ground states IS with 
the same quark content. The errors in the plot have been cal- 
culated by combining the errors listed in Table FVl in quadra- 
ture. For the D mesons the PDG numbers for the neutral 
states, which are more accurately determined, are plotted. 



effects of finite volume are small at large sea quark 
masses, there might be a concern that this is not the 
case for the ensemble with the lightest sea quark where 
MttL « 2.9/m. However, finite volume effects have been 
studied for heavy mesons in [6l| and for the ground state 
D meson are found to be smaller than 1% for the pion 
mass and lattice volume used here. Overall the follow- 
ing picture emerges: Given the limitations of the current 
analysis the = ^ and = |^ multiplets are in 
reasonable agreement with experiment, while the large 
masses of the heavy-strange = doublet resulting 
from our calculation are hard to explain by these limita- 
tions. In particular, we consider it unlikely that a contin- 
uum extrapolation will change these results qualitatively. 

In the case of the heavy-strange jP = i states, the 
DK and D*K thresholds may play an important role. 
Figure fT6l shows the lowest energy levels for the D*q (left 
panel) and Dgi (right panel) ground states. The physical 
DK and D*K thresholds are indicated with a red plus. 
The black plus signs show the central value of these scat- 
tering states on the ensembles with the three lightest sea 
quark masses. While the D mesons have been determined 
in this paper, the kaon mass has been taken from j23| . 
The leading discrepancy between the scattering level at 
the physical point and the scattering level on the lattice 
results from the unphysically heavy strange quark mass 
used in our simulations. An important observation is 
that the overlap of regular quark- antiquark interpolating 
fields on multihadron states (which has to be present in 
dynamical simulations) is most likely small and current 
state of the art spectrum calculations [62], [6^ do not see 
most of the possible scattering states with a simple qq ba- 
sis. (For states in a P-wave, the contribution from scat- 
tering states is expected to be volume suppressed. How- 
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Ds2(2573) -T5 


536.7 ±9.2 ±7.5 


496.3 ± 1.0 


D' -IS 


654.4 ±26.7 ±9.2 


- 


D*' -IS 


726.4 ±20.8 ±10.2 


632.7(±9)(-6) 


2S-TSDs 


708.4 ± 19.9 ±9.9 


- 


IS hyperfine Ds 


133.1 ± 1.0 ± 1.9 


143.8 ±0.4 


2S hyperfine Ds 


72 ± 24 ± 1 


- 


Do (2400) - IS 


266.9 ± 17.3 ±3.7 


347 ± 29 


Di(2420) - IS 


399.1 ± 13.5 ± 5.6 


451.6 ± 0.6 


Di(2430) - IS 


525.2 ± 19.4 ±7.4 


456 ± 40 


L»2(2460) - IS 


577.1 ±20.3 ±8.1 


491.4 ± 1.0 


IS hyperfine D 


130.8 ±3.2 ±1.8 


140.65 ±0.1 



TABLE V: Mass differences in the D and Ds spectrum in 
MeV compared to experimental values (calculated from [J]). 
Bars denote spin-averaged values as discussed in Section HVl 
For the results of this paper, the first error denotes the sta- 
tistical error and the second error denotes the error from set- 
ting the lattice scale. In addition there is a non-negligible 
error from the uncertainty in the determination of Kc for all 
spin-dependent quantities. For the hyperfine-splitting in the 
heavy-strange system, this uncertainty is estimated from the 
tuning run to be 3.2MeV. For a better comparison with the 
heavy-strange mesons, the hyperfine-splitting for the charged 
D mesons is used. We stress again that the gauge ensem- 
bles at our disposal do not allow for a continuum and infinite 
volume extrapolation. 



ever, for many channels even the expected S-wave levels 
are absent.) Comparing the eigenvectors at our light- 
est and third-lightest pion mass suggests that we observe 
the same state on all ensembles and that the expected 
S-wave scattering level is also absent from our data. The 
obvious solution is to include multihadron states in the 
basis for the variational method, which is a challenging 
task left for future calculations. Including multihadron 
states in the basis would also enable one to treat these 
mesons correctly as resonances. The mass and width 
of hadronic resonances below the inelastic threshold can 
be determined within Liischer's finite volume framework 
[g^ [65j . First attempts to include scattering states for 
the calculation of t he p resonance mass and width have 
been published in |66l - l69| and similar studies are cur- 
rently underway for other light-quark resonances. For 
charmonium excitations, multiparticle states have been 
included in [tO] ■ For heavy-li ght mesons first steps in this 
direction have been made in j7l| . 

Finally, there is also the possibility that the states 
in question are indeed not conventional quark-antiquark 
states, but rather tetraquark states or other states with 
which our basis of interpolating fields has poor overlap. 
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FIG. 16: Measured energy levels for the D*q (left panel) and 
Dal (right panel) ground states (black circles) compared to 
experimental states (magenta circles). All masses are plotted 
with respect to the spin-averaged Ds ground state. The plus 
signs denote the DK and D*K scattering levels on the lattice 
(black) and in nature (red). At our lowest pion mass the arti- 
ficially heavy scattering states are very close to the measured 
ground state energy. 



For the case of the D*q as a tctraquark, this has been in- 
vestigated in lattice QCD recently [t^ and no low-lying 
tetraquark states have been identified. 
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Appendix A: Tables 

1. Interpolating fields 

Table IVII collects the interpolators used for both char- 
monium and heavy-light mesons. For the details of the 
source construction, please refer to Section fll Al 



2. Results 

In this appendix, the combination of interpolating 
fields used for the fits, the respective fit intervals 
and results for the energy levels are tabu- 



mm 1 ^max 



[t 

lated in Tables EU IVIlIl and El The ensembles are 
labeled by their number in Table HI We also provide the 
resulting x^/d.o.f. of the correlated fits and the reference 
time slice whenever a matrix of interpolators was used. 
The states are labeled by their quantum numbers. Ex- 
cited states are listed right below the respective ground 
states. For a full list of interpolating fields please refer 
to Section lAll 
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Lattice irrep 


Quantum numbers in irrep 


Interpolator label 


Operator 


Al- 


0-,4-, ... 


A 








B 








C 








D 




A1+ 


0+,4+, ... 


A 


q^q's 






B 


qsiiq'n. 






C 


qs-yniq'D, 


rr 


1-, 3-, 4", ... 


A 


qsjiq's 






B 


qs^t"fiq's 






C 


qsq'oi 






D 


qseijk^j^sq'o^ 






E 


qsltq'o. 






F 


qstijk'yt'YjToq'D,- 




l+,3+,4+, ... 


A 


qsiiToq'a 






B 


qs^iikjiquk 






C 


qstijklfyjq'D^ 






D 


qsjtjz'ysq's 






E 


qsisq'o. 






F 


7*7590, 


T2 


2~, 3-, 4-, ... 


A 


qshjkhjToq'ok 






B 


qs\(-ijk\iti]i5q'Di, 


T+ 


2+, 3+,4+, ... 


A 


qs\eijk\'yjq'D^ 






B 


gs|eijfe|7t7jtoi. 


E- 


2-, 4-, ... 


A 


qsQijk-yrm'Dk 






B 


qsQijkjt'yjToq'D^ 


E+ 


2+,4+, ... 


A 


qsQijk"fjq'D^ 






B 


qsQijkltljq'o^ 
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reduced lattice symmetry implies an infinite number of continuum spins in each irreducible representation of the octahedral 
group. For operators, repeated roman indices indicate summation. The quantity 74 denotes the Dirac matrix for the time 
direction. The quarks, which are in general not of the same flavor, are referred to by q and q . The subscripts s and Di denote 
Gaussian smeared and derivative quark sources respectively. For charmonium, only a subset of operators has been used. Please 
refer to Section Hill for some cautionary remarks regarding Charmonium and charge conjugation. 
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State 


Interpolators 


Ensemble 


tmin 




Energy 


xVd.o.f. 


to 


^ ground state 


A,B,C 


1 


7 


26 


1.22924(23) 


1.22 


4 






2 


7 


27 


1.22098(21) 


2.03 


4 






3 


7 


27 


1.21369(20) 


0.75 


4 






4 


7 


27 


1.21017(24) 


0.71 


4 






5 


7 


27 


1.21268(15) 


0.71 


4 






6 


7 


27 


1.20765(21) 


1.82 


4 


0""*" first excitation 


A,B,C 


1 


7 


11 


1.5595(142) 


1.21 


4 






2 


7 


15 


1.5341(136) 


0.53 


4 






3 


7 


14 


1.5395(97) 


0.48 


4 






4 


7 


12 


1.5247(122) 


2.12 


4 






5 


6 


14 


1.5385(71) 


1.78 


4 






6 


6 


13 


1.5151(82) 


0.59 


4 


1 ground state 


A,B,C,D 


1 


7 


26 


1.27877(37) 


1.01 


4 






2 


7 


27 


1.26860(33) 


2.08 


4 






3 


7 


27 


1.25988(31) 


1.39 


4 






4 


7 


27 


1.25499(34) 


0.99 


4 






5 


7 


27 


1.25904(22) 


0.93 


4 






6 


7 


27 


1.25258(35) 


1.45 


4 


1^^ first excitation 


A,B,C,D 


1 


7 


15 


1.5876(110) 


0.85 


4 






2 


7 


15 


1.5824(110) 


0.83 


4 






3 


7 


13 


1.5535(95) 


0.09 


4 






4 


7 


12 


1.5467(99) 


3.16 


4 






5 


7 


13 


1.5618(71) 


2.35 


4 






6 


7 


13 


1.5373(86) 


0.83 


4 


0"*"^ ground state 


A,B,C 


1 


7 


24 


1.4597(29) 


0.46 


3 






2 


7 


22 


1.4468(27) 


1.83 


3 






3 


7 


23 


1.4256(29) 


0.57 


3 






4 


7 


23 


1.4223(29) 


0.62 


3 






5 


7 


23 


1.4289(18) 


0.33 


3 






6 


7 


24 


1.4165(28) 


0.70 


3 


ground state 


A 


1 


7 


24 


1.4937(30) 


0.73 


- 






2 


7 


25 


1.4756(30) 


1.41 


- 






3 


7 


24 


1.4584(27) 


0.72 


- 






4 


7 


25 


1.4538(26) 


0.75 


- 






5 


7 


24 


1.4639(22) 


0.69 


- 






6 


7 


26 


1.4475(27) 


0.77 


- 


2"*"^ ground state 


A,B in T2 irrep 


1 


7 


22 


1.5095(34) 


0.97 


3 






2 


7 


22 


1.4961(33) 


0.74 


3 






3 


7 


21 


1.4808(29) 


0.86 


3 






4 


7 


26 


1.4735(27) 


0.96 


3 






5 


7 


24 


1.4856(22) 


0.91 


3 






6 


7 


24 


1.4661(26) 


0.88 


3 


1+^ ground state 


D 


1 


7 


21 


1.4932(31) 


0.80 


- 






2 


7 


22 


1.4826(39) 


0.98 


- 






3 


7 


24 


1.4670(31) 


0.39 








4 


7 


21 


1.4627(34) 


1.19 








5 


7 


24 


1.4711(27) 


0.97 








6 


7 


25 


1.4571(33) 


1.38 





TABLE VII: Tabulated fit results for charmonium. The states are labeled by their quantum numbers . The interpolator 
labels and associated structures can be found in Appendix I A II The ensembles are named as in Table [T] The energy levels and 
the associated x^/d.o.f. are from fully correlated two parameter fits to the eigenvalues of the generalized eigenvalue problem 
(GEVP) or to single diagonal correlators. Where applicable, the reference time to for the GEVP is also tabulated. 



17 



State 


Interpolators 


Ensemble 


train 


tmax 


Energy 


xVd.o.f. 


to 


0^ ground state 


A,B,C 


1 


6 


26 


0.86273(33) 


1.42 


3 






2 


6 


26 


0.85158(37) 


1.04 


3 






3 


6 


24 


0.84000(36) 


1.16 


3 






4 


7 


25 


0.82848(40) 


1.37 


3 






5 


6 


24 


0.83929(26) 


1.22 


3 






6 


7 


25 


0.83149(34) 


1.80 


3 


0^ first excitation 


A,B,C 


1 


6 


11 


1.2423(175) 


1.69 


3 






2 


6 


11 


1.2361(205) 


1.74 


3 






3 


6 


11 


1.2406(157) 


1.19 


3 






4 


6 


12 


1.2418(168) 


1.65 


3 






5 


6 


12 


1.2294(117) 


1.32 


3 






6 


6 


13 


1.1782(123) 


1.24 


3 


1~ ground state 


A,B,C,D 


1 


6 


28 


0.93297(59) 


0.73 


3 






2 


6 


26 


0.91864(63) 


1.51 


3 






3 


6 


25 


0.90429(60) 


1.57 


3 






4 


7 


27 


0.89015(69) 


1.20 


3 






5 


6 


24 


0.90429(43) 


0.68 


3 






6 


7 


30 


0.89268(59) 


1.88 


3 


1~ first excitation 


A,B,C,D 


1 


6 


13 


1.2772(135) 


0.63 


3 






2 


6 


10 


1.2732(150) 


0.91 


3 






3 


6 


12 


1.2600(126) 


0.95 


4 






4 


6 


12 


1.2794(117) 


2.07 


3 






5 


6 


12 


1.2545(88) 


1.03 


3 






6 


6 


13 


1.2113(96) 


0.56 


4 


0^ ground state 


A,B,C 


1 


6 


21 


1.1057(35) 


0.98 


3 






2 


6 


21 


1.0784(40) 


0.59 


3 






3 


6 


20 


1.0509(43) 


0.92 


3 






4 


7 


13 


1.0399(44) 


1.28 


3 






5 


6 


24 


1.0570(25) 


0.68 


3 






6 


6 


20 


1.0342(37) 


1.54 


3 


1+ ground state 


A,C,D,F 


1 


7 


13 


1.1574(35) 


0.12 


4 






2 


7 


15 


1.1319(35) 


1.22 


4 






3 


7 


15 


1.1073(34) 


1.01 


4 






4 


7 


13 


1.0907(33) 


2.63 


4 






5 


7 


15 


1.1130(23) 


1.04 


4 






6 


7 


14 


1.0888(31) 


0.79 


4 


1'*' first excitation 


A,C,D,F 


1 


7 


13 


1.1654(54) 


0.82 


4 






2 


7 


15 


1.1458(55) 


1.02 


4 






3 


7 


15 


1.1273(49) 


0.27 


4 






4 


7 


13 


1.1155(51) 


2.92 


4 






5 


7 


15 


1.1320(35) 


1.22 


4 






6 


7 


14 


1.1047(44) 


1.01 


4 


2+ ground state 


A,B in T2 irrep 


1 


7 


20 


1.1810(52) 


0.76 


3 






2 


7 


17 


1.1662(60) 


1.11 


3 






3 


7 


18 


1.1491(52) 


0.72 


3 






4 


7 


18 


1.1353(51) 


1.47 


3 






5 


7 


19 


1.1530(36) 


1.41 


3 






6 


7 


20 


1.1241(44) 


1.29 


3 



TABLE VIII: Tabulated fit results for the Ds mesons. The states are labeled by their quantum numbers . The interpolator 
labels and associated structures can be found in Appendix lA II The ensembles are named as in Table ID The energy levels and 
the associated x^/d.o.f. are from fully correlated two parameter fits to the eigenvalues of the generalized eigenvalue problem 
(GEVP) or to single diagonal correlators. Where applicable, the reference time to for the GEVP is also tabulated. 
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State 


Interpolators 


Ensemble 




tmax 


Energy 


xVd.o.f. 


to 


0~ ground state 


A,B,C 


1 


6 


26 


0.84022(37) 


1.52 


3 






3 


7 


23 


0.79580(61) 


0.79 


3 






5 


7 


26 


0.78798(82) 


1.35 


3 






6 


7 


24 


0.77646(119) 


1.23 


3 


1~ ground state 


A,B,C,D 


1 


5 


26 


0.91237(60) 


0.89 


3 






3 


7 


26 


0.86327(99) 


0.92 


3 






5 


7 


22 


0.85776(122) 


1.41 


3 






6 


7 


24 


0.83656(189) 


1.10 


3 


0^ ground state 


A,B,C 


1 


6 


15 


1.0911(50) 


0.28 


3 






3 


6 


11 


0.9956(80) 


0.38 


3 






5 


6 


17 


1.0150(66) 


1.04 


3 






6 


6 


18 


0.9442(79) 


0.68 


3 


1^ ground state 


A,D 


1 


7 


14 


1.1394(47) 


1.09 


4 






3 


7 


15 


1.0631(48) 


0.64 


4 






5 


7 


12 


1.0601(58) 


2.07 


4 






6 


7 


15 


1.0050(63) 


0.84 


4 


l"*" first excitation 


A,D 


1 


7 


14 


1.1538(65) 


1.00 


4 






3 


7 


15 


1.0874(69) 


0.77 


4 






5 


7 


12 


1.0883(78) 


0.65 


4 






6 


7 


15 


1.0629(93) 


0.79 


4 


2+ ground state 


A,B in T2 irrep 


1 


7 


18 


1.1869(60) 


0.85 


3 






6 


7 


18 


1.1153(74) 


1.48 


3 






3 


7 


15 


1.1166(80) 


1.02 


3 






5 


7 


18 


1.0868(95) 


1.10 


3 



TABLE IX: Tabulated fit results for the D mesons. The states are labeled by their quantum numbers . The interpolator 
labels and associated structures can be found in Appendix I A II The ensembles are named as in Table [T] The energy levels and 
the associated x^/d.o.f. are from fully correlated two parameter fits to the eigenvalues of the generalized eigenvalue problem 
(GEVP) or to single diagonal correlators. Where applicable, the reference time to for the GEVP is also tabulated. 
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